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1 INTRODUCTION 



ABSTRACT 

Comparison of the observed evolution of the Ly-alpha transmitted flux in the spectra 
of four highest redshift quasars discovered by SLOAN survey with the theoretical 
prediction for this evolution based on the state-of-the-art numerical simulations of 
cosmological reionization already allows one to constrain the redshift of reionization 
to zrei = 6.2 ± 0.1s ±0.2r, where systematic and random errors are given respectively. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - galax- 
ies: formation - galaxies: intergalactic medium 



Table 1. Simulation Parameters 



In a recent paper Becker et al. (2001) analysed spectra of the 

■ four highest redshift quasars discovered by SLOAN survey. 
The highest redshift (z — 6.28) quasar showed almost no 
transmitted flux just blueward of the quasa r Ly- alpha emis- 

■ sion line, and on this bases Becker et al. (2001) concluded 
that reionization took place at a redshi ft z ~ 6 [a similar 



claim was also made by Djorgovski et al. (2001) based on an 
extrapolation from a lower redshift observations]. 

However, one should be cautious before drawing such a 
conclusion. Indeed, the observed decrease in the mean trans- 
mitted flux at z ~ 6 might simply indicate a decrease in the 
mean ionizing intensity rather then a real reionization of the 
universe. Thus, without an understanding of the evolution 
of the universe around the reionization epoch, the Ly-alpha 
absorption data cannot be used to constrain the epoch of 
reionization [unles s damping wings are observed in the ab- 
sorption profiles, ( Miralda-Escude 1996 )]. 

Fortunately, our theoretical understanding of the pro- 
cess of reionization, in part based on numerical simulations, 
is solid enough so that the evolution of the ionizing intensity 
at these redshifts can be predicted with a reasonable con- 
fidence level. Combining simulations with the observational 
data indeed allows one to come up with a meaningful value 
for the redshift of reionization. 



2 SIMULATIONS 

A set of three simulations have been performed with the 
SLH code and are similar to the simulations reported in 
(Gnedin 2000a). The main difference with previous simula- 



tions is that a newly developed and highly accurate Opti- 
cally Thin Eddington Variable Tensor (O TVET) approx- 



imat ion for modeling radiative transfer (Gnedin & Abel 
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0.15 


2.5 x 10~ 6 



ciently accurate (subject to the usual limitations of numer- 
ical convergence and phenomenological description of star 
formation) to be used meaningfully in comparing with the 
observational data. 

Parameters of the three simulations are given in Table |l[ 
All three simulations included 64 3 dark matter particles, an 
equal number of baryonic cells on a quasi-Lagrangian mov- 
ing mesh, and about 100,000 stellar particles that formed 
continuously during the simulation. The box size was fixed 
at 2/i -1 comoving Mpc, and the nominal spatial resolution 
of the simulation was fixed at l/i -1 comoving kpc, with the 
real resolution being a factor of two worse. 

In all cases a flat cosmology was assumed, with £Ia,o = 
1 — n m> o, and COBE normalization was adopted. Notice 
that because of that, a small change in the slope of the pri- 
mordial power-law spectrum n makes a significant effect on 
the amount of the small-scale power due to a large leverage 
arm from COBE scales to the tens-of-kpc scales which are 
important for reionization. 

Star formation is incorporated in the simulations using 
a phenomenological Schmidt law, which introduces two free 
parameters: the sta r form ation efficiency csf [as defined by 
eq. (6) of Gnedin (2000b)] and the ionizing radiation effi- 



2001) is used instead of a crude Local Optical Depth ap- 
proximation. The new simulations therefore should be suffi- 



ciency euv (defined as the energy in ionizing photons per 
unit of the rest energy of stellar particles). These two pa- 
rameters are listed in the last two columns of Table [I] for 
reference purposes. 

The star formation efficiency esf is chosen so as to 
normalize the global star formation rate in the simula- 
t ion a t z — 4 to the observed value from Steidel et al. 
(h_999(), whereas the ultraviolet radiation efficiency euv is 
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Figure 2. Mean transmitted flux as a function of redshift for 
simul ations A-C. Points show the observational data from Becker 
et al. (2001). The last two data points (open circles) are obtained 



from only one quasar, and the vertical error due to cosmic vari- 
ance cannot be estimated. Arrows mark the redshift of reioniza- 
tion for each simulation. 



Figure 1. Mean free path to ionizing radiation (bottom) and 
its time derivative (top) as a function of redshift for the three 
simulations from Table 



only weakly constrained by the (highly uncertain) mean 
photoionization rate at z ~ 4. The redshift of reionization 
strongly depends on euv and is not in fact predicted in a 
simulation, but can be changed over a reasonable range de- 
pending on the assumed value of euv • I will use th is fre edom 
to fit the observational data from Becker et al. (2001) and 
instead to constrain the redshift of reionization from the 
observational data. 



A "REDSHIFT OF REIONIZATION": 
IS IT? 



WHAT 



Before one can attempt to measure "the redshift of reion- 
ization" , we better be sure that such a quantity can be de- 
fined. The whole process of reionization is quite extended 
(Az ~ 5 — 10), and even the fast process of percolation of 
ionized bubbles occurs over a sizable redshift interval Az ~ 1 



(Gnedin 2000a). However, one can still define the specific 



value of the redshift of reionization as the moment which 
corresponds to the peak rate of increase of the mean free 
path to ionizing radiation. As can be seen from Figure |jj 
the time derivative of the mean free path has a well defined 
peak, which I use throughout this paper as "the redshift of 
reionization" zrei- 



4 RESULTS 

Figure ^ shows the mean transmitted flux as a function of 
redshift for the three simulations described above. In each 
case the arrow shows the redshift of reionizatio n, and points 
with error-bars are taken from Becker et al. (2001). Each 
point represents an average of up to four quasars, except for 
the last two points that are derived from a single quasar and 



thus have no vertical error-bars. It is important to empha- 
size that the measurement of the mean transmitted flux for 
each quasar over an interval Az ~ 0.15 is quite accurate, 
with the intrinsic error of only 0.003 — 0.005 in (F), but 
the variation between different lines of sight (the so-called 
"cosmic variance") is much greater, and it is this variation 
that dominates the vertical error-bar. Thus, one (or better 
two) more z = 6.3 quasars are required in order to place the 
vertical error-bars on the open circles. However, because the 
sharp drop in the mean transmitted flux at z ~ 6 is marked 
by two data points, it is more reliable than simply one point 
from one quasar. 

Simulations in Fig. |i| differ from the data points in two 
ways: both the redshift evolution in a simulation and the 
photoionization rate after reionization (the amplitude of the 
curve at low z) are offset relative to the data. As I mentioned 
above, the free parameter euv can be used to adjust the 
simulation to fit the observational data. However, there is no 
guarantee that with one parameter I can adjust two offsets 
at the same time for a given cosmological model. This fact 
is extremely important because it allows one to actually put 
constraints on the cosmological model per se, and I elaborate 
on this opportunity in the conclusions, but here I am going 
to ignore this fact and adjust two offsets independently - 
by sliding the curve both vertically and horizontally - to fit 
the observational data. Because the three curves from three 
simulations have similar shapes, every simulation can thus 
be made to fit the data - and, again, in reality, only a narrow 
range of cosmological models will succeed in doing so. 

The reason for doing so is to obtain a constraint on 
the redshift of reionization which does not depend on a 
(weakly constrained) cosmological model. In addition, the 
simulations presented here are rather small and numeri- 
cal errors due to incomplete convergence are substantial 
( |Gnedin 2000al ). Simulat ions with larger box sizes typically 
have lower mean transmitted flux after reionization than 
small box simulation - which implies that amplitudes of 
three curves are not sufficiently accurate in Fig.]2[ The ar- 
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larger than a factor of 3-5. One or two more z = 6.3 quasars 
are required to confirm or refute this assumption. 

A strong sensitivity of the redshift of reionization to 
the amount of small-scale power offers a unique opportunity 
to place constraints on the slope of the primordial power 
spectrum - unattainable even with the SLOAN data on the 
power spectrum at megaparsec scales. For example, runs 
B and C in Fig. ^ have similar star formation rates and 
redshifts of reionization, but differ by a factor of 3 in the 
mean transmitted flux at z — 5 — 5.5 only because the slope 
of the primordial power spectrum n differs by mere 0.02 
in the two models. Sufficiently large simulations that have 
numerical effects under control are currently feasible and 
will eventually provide a tight constraint on the amount of 
small-scale power. 



Figure 3. Mean transmitted flux as a function of redshift for 
sim ulatio ns A-C, fitted to the observational data from Becker et 
al. ( |200l[ ) by shifting curves both in vertical and in horizontal 
direction. Shaded regions illustrate the uncertainty due to obser- 
vational errors; this uncertainty is also given by the error-bars 
drawn on the arrowheads. 



bitrary vertical shift of the curves can thus be considered as 
a "marginalization" over the simulation box size. 

Figure ^ shows the fit to the observational data for each 
of the three simulations. The shaded region also gives an 
uncertainty of this fit, which can be considered a random 
uncertainty due to observational errors (mostly cosmic vari- 
ance). Because the last two points have no vertical error- 
bars, they are only partly used: they only constrain models 
in the horizontal direction, and their redshift error-bars un- 
certainty are somewhat arbitrarily increased by a factor of 
2 (to account for the possibility that they can be moved up 
and down). The fact that different cosmological models have 
slightly different redshifts of reionization when made to fit 
the observed evolution of the mean transmitted flux illus- 
trates the systematic uncertainty due to unknown cosmo- 
logical parameters. With these two uncertainties included, 
I can derive a value for the redshift of reionization in the 
following form: 



zrei = 6.2 ± 0.1s ±0.2r, 



(1) 



which is the sole result of this paper. 

If cosmic variance on the last two points can be esti- 
mated, and if it is comparable to the vertical error-bar on 
the z = 5.7 point, then the random error gets reduced from 
0.2 to 0.1. 
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5 CONCLUSIONS 

SLOAN observations of z ~ 6 quasars push the frontier of 
the observable universe right into the epoch of reionization. 
Combined with the most advanced simulations of cosmo- 
logical reionization to date, the observational data yield a 
rather precise measurement of the redshift of reionization. 

This measurement hinges on the assumption that the 
z — 6.28 quasar (SDSSp 1030+0524) probes an average re- 
gion of the universe, i.e. that the cosmic variance in the 
measurement of the mean transmitted flux at z — 6.1 is not 
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